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A high resolution UV absorption spectrum of supernova 
ejecta in SN1006 
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ABSTRACT 

We report a high resolution, far- ultraviolet, STIS E140M spectrum of the strong, broad 
Si II, III, and IV features produced by the ejecta of SN1006 seen in absorption against 
the background Schweizer-Middleditch star. The spectrum confirms the extreme 
sharpness of the red edge of the redshifted Si II 1260 A feature, supporting the idea 
that this edge represents the location of the reverse shock moving into the freely 
expanding ejecta. The expansion velocity of ejecta at the reverse shock is measured to 
be 7026 ± 3 (relative) ± 10 (absolute) km s" 1 . If the shock model is correct, then the 
expansion velocity should be decreasing at the observable rate of 2.7 ± 0.1 km s _1 per 
year. The pre-shock velocity, post-shock velocity, and post-shock velocity dispersion 
are all measured from the Si II 1260 A feature, and consistency of these velocities with 
the shock jump conditions implies that there is little or no electron heating in this 
fast (2680 km s" 1 ) Si-rich shock. 
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1 INTRODUCTION 

The supernova o f 1006 AD (SN1006) was t he brightest in 
recorded history dStephenson fc Greenll20o3l . It is generally 
held that SN1006 was a probable Type la, the most straight- 
forward evidence for this being its location high above the 
plane of the galaxy, ~ 500 pc at the 2 . 18 ±0 .08 kpc distance 
measured bv lWinkler. Gupta fc Lonel J2003f) . The identifica- 
tion of SN1006 as Type la is consistent with a broad range 
of observational evidence and theoretical expectation. The 
historical supernova was exceptionally bright, and remained 
visible over an extended (3 year) period. There is no sign of 
a pulsar or other central compact object. 

The most likely cause of a Type la supernova is thought 
to be the thermonuclear explosion of a white dwarf which ac- 
cretes from a companion, probably a giant or main sequence 
star, but possibly a wh i te dwarf l|Hillebrandt fc Niemeverl 
l2000l:lHoflich et all20o3lHoflicrl200^ . Searches have so far 
failed to uncover any candidate star that may have been th e 
erstwhile companion in SN1006 jRuiz-Laouente et all20o3i . 

The remnant of SN1 006 has been ext ensively 
observed in the radio (Reynolds fc Gilmord 
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IDver. Reynolds fc Borkowskll 12004 IVinkl [ 20051 
Jiaienic^^]tLl~|200fi) a nd o '-ravs iTmhnoTietHil] ^98l . 
2o7nT r TAhliro'riian"eta^ 120051) . The observations show a 
limb-brightened, roughly circular blast wave expanding into 
an interstellar medium of moderate density, ~ 0.1 cm" 3 . 
The x-ray image and spectrum reveal two components 
of emission: limb-concentrated non-thermal synchrotron 
radiation from relativistic electrons accelerated at the 
blast wave, and more generally distributed lumpy thermal 
emission from a hot plasma enhanced in heavy elements 
jLong et all l200l Fig. 2). The thermal x-ray emission is 
consistent with what might b e expected from t he ejecta of a 
Type la supernova explosion iDver et all200ll . Table 5) , but 
the inferences about masses and elemental abundances are 
model dependent, and alternative inte rpretations cannot be 
excluded (but see lBadenes et all | 2u 03 for a more optimistic 
view) . 

One of the more direct observational evidences on 
the nature of SN1006 comes from absorption spec- 
tro scopy of a background OB su bdwarf star discovered 
by ISchweizer fc Middleditchl (Il980l hereafter the SM star) 
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Table 1. STIS Observing log 

Object RA & Dec (J2000) Date Orbits Exposure (s) 

PG 0839+399 08 h 43 m 12.70 s +39°44'49'.'9 21 Oct 1998 1 2,494 

SM star 15 h 02 m 53.18 s -41°59'17'.'6 24-25 Jul 1999 8 2 x (2,527 + 9,066) = 23,186 



just 2.8' south of the projected centre of the 15' radius 
remnant of SN1006. Ultraviolet spectroscopy of the SM 
star with the In ternational Ultraviolet Explorer (IUE) by 
revealed broad absorption features, which 



Wu et alJ 11983ft attributed to Si II, 
Fe II in the ejecta of SN1006. The 



III, and 



IV, and 
of these 



presence 

featu res was confirmed in s ubsequent observations with 
IUE iFesen fe Hamiltonl H988). and at higher resolution and 
signal-to- noise with the Faint Object Spe ctrograph (FOS) 
on t he Hubble Space Telescop e (HS T) l|Wu et alJ |l993t 
Il997ft . iBlair. Long fe Raymond! lll996ft used the Hopkins 
Ultraviolet Telescope (HUT) to observe the spectrum short- 
ward of Lya, finding marginal evidence for the presence of 
absorption by Fe III 1123 A. Optically, the SM star shows 
no ev idence for absorption by SN1006 ejecta feurleigh et alJ 
12000ft . Recently. IWinkler et aD (12005ft have reported obser- 
vations of four more UV sources behind SN1006, finding 
clear evidence of absorption by SN1006 ejecta in the two 
sources nearer to the projected centre. 

In an exte nsive analysis of the H ST FOS spectrum 
of the SM star, lHamilton et al.l (|l997j, hereafter Paper 1) 
argued that the unusual shape of the broad Si II 1260 A 
feature, the strongest absorption feature in the 1150-2000 A 
ultraviolet spectrum of the SM star, arises because the 
feature consists of both unshocked and shocked components. 
This surprising conclusion was driven by several observa- 
tional facts. (1) The Si II 1260 A feature had a sharp red 
edge, at 7070 ± 50 km s" 1 , which Paper 1 identified as the 
position of the reverse shock front in the freely expanding 
ejecta of SN1006. This suggested that the reverse shock in 
SN1006 is currently moving into ejecta containing Si II. (2) 
The depth of absorption of the Si II 1260 A feature implied 
that the collisional ionization time of shocked Si II is of 
the order of the 1000 yr age of the remnant. This implied 
that there should also be Si II in the shocked ejecta. (3) 
The blue edge of the Si II 1260 A feature fitted nicely to a 
broad Gaussian profile, as expected for shocked Si II, and 
the same Gaussian profile accounted naturally for the slight 
tail of absorption observed to high velocities > 7000 km s _1 . 
(4) The velocities of unshocked and shocked gas measured 
from the observed Si II 1260 A line profile were consistent 
with the shock jump conditions. 

In this paper we report Hubble Space Telescope Imaging 
Spectrometer (STIS) observations of the SM star behind 
SN1006, with a spectral resolution of 5 km s _1 FWHM, 
~ 60 times the resolution of the earlier FOS spectrum 
jWu et al.lll997ft . A primary motivation for obtaining the 
high resolution STIS spectrum was to test the proposal of 
Paper 1 that the Si II 1260 A feature indeed reveals ejecta in 
SN1006 undergoing a reverse shock. A key prediction of the 
shock model is that the red edge of the Si II 1260 A feature, 
representing the position of the reverse shock, should be 
extremely sharp, since a shock should decelerate gas 'instan- 
taneously' from the free expansion velocity to some lower 
bulk velocity. High resolution spectroscopy can also test for 



possible problems, such as saturation, or contamination by 
narrow lines. 

One of the main sources of unquantified system- 
atic uncertainty in interpreting the absorption features in 
the FOS spectrum was the possibility of contamination 
by stellar lines. We therefore obtained and report here 
STIS observations of a comp arison star PG 0839+399 
jGreen. Schmidt fe Liebertlll986ft identified by IBlair et alJ 
( 1996) as having not only a similar temperature and gravity, 
but also similar photospheric abundances, and low extinc- 
tion. 

The STIS spectra reported here do not cover the Fe II 
2383 A and 26 00 A absorption features observed in the SM 
star with I UE dWu et al]ll983tlFesen fe Hamiltorll98Sft and 
the FOS <Wu et al]ll99l Il997ft . and more recently in a 
second background source, the quas ar QSO-1503-4155, wi th 
low-resolution STIS observations bv I Winkler et all ((2005) . 



2 OBSERVATIONS 

Far-ultraviolet echelle spectra of the SM star and of the 
comparison star PG 0839+399 were obtained in 1998 and 
1999 with STIS using the 2" x 2" slit, the E140M grating, 
and the FUV-MAMA detector. The observing log is given in 
Tabled The spectra cover 1150-1700 A in wavelength bins 
AA ps A/92000, or 3.3 km s _1 , wide. Over the relevant range 
of the spectra from 1220-1540 A, the spectral resolution is 
1.4 pixels, or 4.6 km s _1 , FWHM. The line spread function is 
non-Gaussian, with a narrow core and broader wings. Over 
1220-1540 A, the signal-to-noise ratio is ~ 8 per pixel and ~ 
15 per pixel for the SM star and PG 0839+399 respectively. 

Figure □ shows the STIS E140M spectra of the SM star 
and of the comparison star PG 0839+399. Each spectrum 
has been convolved with a near-Gaussian with a FWHM of 
320 km s _1 . The near-Gaussian is a function proportional 



to (1 



for | ac | < 1, with n = 4. Such a smoothing 



function mostly shares the nice theoretical properties of 
Gaussian smoothing, but extends over only a finite number 
of pixels. The smoothing is broad so as to bring out the 
broad absorption features in the spectrum of the SM star; 
higher resolution versions of the spectra are shown below 
in There are a few narrow gaps in the spectra longward 
of 1600 A, and the spectra were linearly interpolated across 
these before being smoothed. 

The small scale structure in the various spectra in 
Figure Q is not noise, but real, mostly stellar features. 
This fact can be inferred from the good agreement between 
spectra from the individual exposures of the SM star, 
apparent for example in the high-resolution zoom of the Si II 
1260 A feature shown in Figure [3] below. 

The stellar lines of the comparison star PG 0839+399 
are blueshifted, by —32 km s~\ relative to those of the SM 
star. Narrow interstellar absorption lines in PG 0839+399 
show no discernible shift relative to those in the SM star. 

Figure [5] shows the ratio of the spectrum of the SM 
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Figure 1. STIS E140M spectra of the SM star behind SN1006 and of the comparison star PG 0839+399 (the latter multiplied by a 
factor 0.1 to bring it on to the graph). Both spectra have been smoothed to a resolution of 320 km s — 1 FWHM. The smooth lines are 
fitted continua, quintic polynomials in logF as a function of 1/A. Narrow interstellar lines are marked. 
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Figure 2. Absorption spectrum of the SM star, showing the broad rcdshiftcd absorption features of Si II, III, and IV produced by 
supernova ejecta in SN1006. The absorption spectrum is the ratio of interstellar-line-excised, smoothed, continuum-corrected E140M 
spectra of the SM star and of the comparison star PG 0839+399, the latter shown at top. The smooth line passing through the 
absorption spectrum is a model in which the weaker Si II 1304 A and Si II 1527 A features are assumed to have the same profile as the 
strong Si II 1260 A feature, which is assumed to be optically thin, while the Si III 1206 A and Si IV 1394, 1403 A features are assumed 
to have the same Gaussian profile as that of the shocked component of Si II 1260 A, as described in i|3.3l The bottom panel shows the 
residual absorption spectrum of the SM star after dividing by the model spectrum. 



4 A. J. S. Hamilton et al. 



Table 2. Measured parameters 

Parameter STIS (this paper) FOS (Paper 1) 

Expansion velocity r s /t into reverse shock 7026 ± 3 (relative) ± 10 (absolute) km s" 1 7070 ± 50 km s _1 

Mean velocity of shocked Si II 5160 ± 70 km s" 1 5050 ±60 km s" 1 

Dispersion a of shocked Si II 1160 ± 50 km s" 1 1240 ± 40 km s" 1 

Reverse shock velocity v s 2680 ± 120 km s" 1 2860 ± 100 km s" 1 

Shocked Si II : Si III : Si IV 1 : 0.45 ± 0.02 : 0.32 ± 0.02 1 ± 0.03 : 0.43 ± 0.02 : 0.41 ± 0.02 
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Figure 3. STIS E140M spectra of the redshifted Si II 1260 A feature, showing separately the two individual 9066 second exposures 
of the SM star, and the comparison star PG 0839+399. The spectrum of PG 0839+399 has been shifted by +32 km s" 1 to mesh its 
stellar lines with those of the SM star, and divided by 15 to bring its continuum to approximately same level as the SM star. All spectra 
have been smoothed with a near- Gaussian of FWHM 6.8 km s" 1 . The vertical arrow shows the position of the putative reverse shock 
at 7026 ± 10 km s — 1 . The good agreement between the two individual exposures of the SM star demonstrates that the structure in the 
spectra is real, not noise. Sadly, the spectrum of the comparison star is not as similar in detail to that of the SM star as one might 
have liked. 



star to that of the comparison star PG 0839+399. Before 
their ratio was taken, both spectra were linearly interpolated 
across narrow interstellar resonance lines, the spectrum of 
PG 0839+399 was redshifted by +32 km s" 1 to mesh its 
stellar lines with those of the SM star, and both spectra were 
smoothed to a resolution of 320 km s _1 FWHM, the same 
as in Figure The continuum was then removed by fitting 
unabsorbed parts of the ratio spectrum to a quintic polyno- 
mial in log flux versus inverse wavelength. The continuum 
fitting is intended to remove differences both in the intrinsic 
continua of the spectra and in the extinction. We tried 
co rrecting the contin ua for the extinction alone, measured 
bv lBlair et all Jl99rl Table 2) to be E B -v = 0.119 ± 0.005 
for the SM star and E B -v = 0. 045+0.005 for PG 0839+399 
using the extinction curve of ICardelli. Clayton fc Mathisl 
(1989), but this correction leaves noticeable differences in 
the continua, particularly at shorter wavelengths, so we 
preferred the empirical polynomial fit. 



The ratio spectrum of Figure |3 
features, which, following IWu et al. 


shows broad absorption 
(1983) and subsequent 


authors jFesen & Hamilton! Il98ct 


Wu et al. 


1993 




Il997t 



and IV in the ejecta of SN1006. 

As in the individual spectra shown in FigureQ the small 
scale structure in the ratio spectrum of Figure|5|is not noise, 
but real, mostly stellar features. It is thus apparent that the 



spectral type and elemental abundances of the comparison 
star PG 0839+399 are not identical to those of the SM star. 

The smooth line passing through the absorption spec- 
trum in Figure|5|is a model in which the weaker Si II 1304 A 
and Si II 1527 A features are assumed to have the same 
profile as the strong Si II 1260 A feature, with strengths in 
proportion to their wavelengths times oscillator strengths, as 
expected if the lines are optically thin. The Si III 1206 A and 
Si IV 1394, 1403 A features are assumed to have Gaussian 
profiles with the same width as that measured in H3.ll for 
the shocked component of the Si II 1260 A feature. 



3 RESULTS 

3.1 Si II 1260 A feature 

Figure [3] shows the STIS spectrum of the SM star at high 
resolution in the vicinity of the redshifted Si II 1260 A 
absorption feature. The Figure shows individually the two 
long exposure spectra, 9066 s each, of the SM star, and 
the spectrum of the comparison star PG 0839+399, the last 
shifted by +32 km s _1 to mesh its stellar lines with those 
of the SM star, and divided by 15 to bring its continuum 
to approximately the same level as that of the SM star. 
Each spectrum has been smoothed with a near Gaussian of 
6.8 km s" 1 FWHM, just slightly greater than the 4.6 km s" 1 
intrinsic resolution of the spectrum. 
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The arrow at 7026 km s _1 marks the velocity of the 
sharp red edge of the feature, the putative position of the 
reverse shock moving into the freely expanding Si II ejecta. 

The spectrum would seem to confirm the key prediction 
of the shock hypothesis, that the edge of the feature should 
be extremely sharp. At the steepest point, the flux shown in 
Figure|^]rises over an interval of 7 pixels from an absorption 
fraction of ~ 0.3 at 7016 km s _1 to an absorption fraction 
of ~ 0.8 at 7036 km s _1 . We take the position of the reverse 
shock to be 7026 km s~\ at the centre of the steep interval, 
and this is the value given in Table [5] 

We have compared the shape of the steep interval to 
that expected for a step function smoothed over the line 
spread function of STIS with the E140M grating and 2" x 2" 
slit, and conclude that the edge is consistent with being 
unresolved. 

It is natural to assign the uncertainty of the location of 
the reverse shock to be ±10 km s -1 , since that is the extent 
of the steep interval. However, if at some time in the future a 
second epoch observation is taken at similar resolution and 
signal-to-noise, then by comparing the line profiles near the 
shock, it should be possible to measure the change in the 
position of the shock with significantly higher resolution. 
Given the quality of the present spectrum, we judge that 
it should be possible to measure the change in the position 
of the shock to better than a pixel, about 3 km s _1 . We 
therefore choose to quote in Table |5] two separate errors on 
the 7026 km s _1 position of the shock, a relative eror of 
±3 km s~\ which is an estimate of the uncertainty in the 
relative position of the shock between two observations of 
similar quality at two different times, and an absolute error 
of ±10 km s _1 , which is an estimate of the uncertainty in the 
absolute position of the shock front in the freely expanding 
ejecta. 

Figure|^|also supports a second proposition of the shock 
model, that the Si II 1260 A absorption feature is optically 
thin. If the feature were produced by de nse, optically thick 
kno ts of ejecta, a possib i lity d iscussed bv lFesen et, alJ ljl98Sh 
and lFesen fe Hamilton! il988tl . then the absorption feature 
would be choppier, hitting zero at various wavelengths. The 
deepest part of the absorption trough, over 1280-1290 A, 
remains clear above zero. 

Figure [I] shows the absorption spectrum of SN1006, 
the ratio of the spectrum of the SM star to that of 
the comparison star PG 0839±399, in the vicinity of the 
Si II 1260 A absorption feature. The Figure is essentially 
a zoomed-in version of Figure |5] but at four times the 
spectral resolution. As with Figure|2| both spectra were first 
linearly interpolated across narrow interstellar resonance 
lines, and the spectrum of PG 0839±399 was redshifted 
by ±32 km s" 1 . Both spectra were smoothed with a near 
Gaussian of 80 km s _1 FWHM before their ratio was 
taken, and the ratio was then corrected by the same quintic 
polynomial as in Figure [2] 

The maximum velocity 7026 km s _1 of unshocked Si II 
represents the position of the reverse shock, and as discussed 
in Paper 1, Si II passing through the reverse shock will 
be decelerated suddenly from the free expansion velocity 
of 7026 km s _1 to a lower bulk velocity. The collisional 
ionization time of the shocked Si II is comparable to the age 
of the remnant, so it is expected that Si II should persist 
after being shocked. The shocked Si II should contribute a 
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Figure 4. Absorption spectrum of SN1006 near the Si II 1260 A 
feature, showing the best fit Gaussian profile of shocked Si II, 
and the residual unshocked Si II. The spectrum is a higher- 
resolution zoomed-in version of the spectrum shown in Figure |21 
The spectrum is the ratio of interstellar-line-excised, continuum- 
corrected spectra of the SM star and of the comparision star PG 
0839+399, each smoothed to a resolution of 80 km s" 1 FWHM 
before their ratio was taken. The upper spectrum shows the 
interstellar-line-excised, continuum-corrected comparison stellar 
spectrum of PG 0839+399 at the same resolution, shifted by 
+32 km s _1 to mesh its stellar lines with those of the SM star, 
and offset to separate it from the SN1006 spectrum. 



broad Gaussian profile of absorption. Fitting the blue edge 
of the Si II 1260 A feature to a Gaussian, we find the best fit 
shown in Figure [I] The best fit parameters of the Gaussian 
fit are a mean of 5160 ± 70 km s _1 , and a dispersion of 
1160 ± 50 km s" 1 , as listed in Tabled 

Table H also gives, for reference, values previously 
measured from the FOS data in Paper 1. The FOS velocities 
listed in Table |5] are as given by Paper 1, and have not 
been corrected for the —41 km s _1 shift of the FOS 
spectrum measured in H3.2l below. It is to be noted that the 
uncertainties on the STIS parameters are, except for the 
much more precise measurement of the expansion velocity 
into the reverse shock, not much different from the quoted 
uncertainties on the FOS parameters. As in Paper 1, the 
stated uncertainties are formal 1 standard deviation statis- 
tical errors, and do not include the systematic uncertainty 
associated with the stellar spectrum or with placement of the 
continuum. In the present paper the systematics are under 
better control, thanks to the template stellar spectrum from 
the comparison star PG 0839±399, so that, notwithstanding 
the comparable statistical errors, the STIS values in Tabled 
are more reliable, and are to be preferred. 

The consistency of the measured velocities can be 
checked, as in Paper 1, against the jump conditions for 
a strong shock. The shock jump conditions predict that 
the three-dimensional velocity dispersion 3 1//2 <r of the ions 
should be related to the deceleration Av of the shocked gas 
by energy conservation 
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Av 



(1) 



provided that all the shock energy goes into ions. From the 
measured 1-dimensional dispersion a = 1160 km s _1 given 
in Table |21 the observed 3-dimensional dispersion is 

3 1/2 cr = 3 1/2 x (1160 ± 50 km s _1 ) = 2010 ± 90 km s _1 .(2) 

The observed deceleration is the difference between the mea- 
sured free expansion velocity 7026 km s _1 of the unshocked 
gas at the reverse shock front, and the measured mean 
velocity 5160 km s~ x of the Gaussian velocity distribution 
of shocked gas, which is 



Av = (7026 ± 10) - (5160 ± 70) = 1870 ± 70 km s" 



(3) 



As previously found in Paper 1, these values are in good 
agreement (the observed dispersion exceeds that predicted 
from the deceleration by 1.2 standard deviations), encour- 
aging the view that the shock interpretation is basically 
correct. 

The reverse shock velocity v 3 corresponding to the 
observed dispersion a is 



= (16/3) 1/2 cr = 2680 ± 120 km s" 1 



(4) 



which is listed in Table [5] 

As remarked above, the relation Q is valid only if all 
the shock energy goes into ions, with negligible energy going 
into heating electrons. If some of the shock energy went into 
electrons, then the predicted ion velocity dispersion a would 
be decreased, which would worsen the agreement with the 
shock jump conditions. Thus, if the shock interpretation of 
the Si II 1260 A feature is accepted, then the measured 
velocities constitute direct evidence that little electron 
heating occurs in shocks of this kind. This is consistent 
with other observational evidence that l ittle electron heating 
occurs in the fast shocks in SN1006 jLaming et alj|l996t 
iGhavamian et alJl2002t Ivink et al]l2003UVinHl200sF 7 



3.2 Evolution of the reverse shock 

As the reverse shock propagates into the Si II ejecta, the 
free expansion velocity r s /t at the position r s of the reverse 
shock should decrease with time t as 

d(r s /t) _ 

- t 



dt 



(5) 



Given the measured shock velocity v s , equation @, and the 
known t = 993 yr age of the remnant at the 1999 epoch 
of observation, equation JSJ predicts that the expansion 
velocity of the reverse shock should be decreasing at a rate 
of 



= -2.7 ±0.1 km s -1 yr" 



(6) 



d(r 3 /t) 
dt 

We have attempted to measure the expected change in 
the position of the reverse shock between the FOS and STIS 
data. The time that elapsed between the FOS observa tions, 
taken during 17-20 September 1993 dWu et al.lll997l) . and 
the STIS observations, taken during 24-25 July 1999, is 
2135 days, or 5.85 years, so if the shock model is correct 
then it is expected that the position of the reverse shock 
should have changed by A(r 3 /t) — 5.85 yr x (—2.7 ± 



0.1 km s 

A{r s /t) 



yr 



" 1 ), that is by 
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-15.8 ± 0.7 km s" 



(7) 
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Figure 5. Comparison of the FOS G130H (thick line) to 
the STIS E140M (thin line) absorption spectra of SN1006 near 
the Si II 1260 A feature. The FOS spectrum has been shifted 
by +41 km s — 1 to bring its stellar and interstellar lines into 
wavelength registration with the STIS spectrum, and adjusted 
by a polynomial to cancel a slowly varying difference in the FOS 
and STIS continua. The STIS spectrum has been convolved to 
the resolution of the FOS spectrum, 370 km s _1 FWHM. Both 
spectra have been divided by the smoothed comparison spectrum 
of PG 0839+399, and corrected by the same quintic polynomial 
applied to the ratio spectrum in Figure 1^1 



between the FOS and STIS observations. 

To measure the change in the position of the reverse 
shock as accurately as possible, we carried out a differential 
comparison of the FOS and S TIS s pectra. The FOS spec- 
trum was taken bv lWu et al] il997f) using the 1" diameter 
circular aperture, the G130H grating, and the blue detector. 
The uncertainties in the comparison are dominated by 
uncertainties in the FOS spectrum. We took the error in each 
FOS spectral bin from the photon-counting noise tabulated 
in the HST spectral data file, multiplied by an empirical 
factor of 1.2 suggested by the scatter in the flux difference 
between adjacent spectral bins. Values and la uncertainties 
of parameters given below are deduced from the minimum 
and deviation in the x 2 01 the difference between the FOS 
and STIS spectra. 

To compare the FOS and STIS spectra, we first 
smoothed the STIS spectrum with a near-Gaussian to the 
resolution of the FOS spectrum. Fitting the smoothed STIS 
spectrum to the FOS spectrum yields a best fit resolution 
of 370 ± 10 km s _1 FWHM, somewhat larger than the 
nominal 295 ± 10 km s _1 FWHM of the line spread function 
of the FOS G130H spectrum in the configuration used. 
Comparison of stellar and interstellar lines in unabsorbed 
(by SN1006 ejecta) ranges of the spectrum over 1240-1400 A 
indicate that the FOS spectrum is shifted by —41 + 6 km s _1 
relative to the STIS spectrum. To bring the FOS lines into 
registration with those of STIS, we therefore shifted the FOS 
spectrum by +41 km s -1 . The FOS and STIS continua differ 
somewhat, the FOS continuum varying from ~ 10% high 
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near 1170 A to about ~ 10% low near 1400 A. To cancel this 
variation, we fitted the log of the FOS to STIS flux ratio to 
a 19th order polynomial as a function of wavelength, and 
divided the FOS spectrum by the fitted ratio. Despite the 
high order, the correcting polynomial varies moderately and 
smoothly with wavelength. 

Figure[S]compares the resulting adjusted FOS spectrum 
to the smoothed STIS spectrum in the vicinity of the Si II 
1260 A feature. 

Having made these preliminary adjustments, we mea- 
sured the change in the velocity of the shock front between 
the FOS and STIS spectra from the shift in wavelength of 
the sharp edge at 1289.5-1291 A. The measured shift of the 
edge of the STIS spectrum relative to the FOS spectrum is 

A(r s /t) = -35 ±25 km s" 1 (8) 

in which approximately half the error comes from uncer- 
tainty in the wavelength of the edge, and the other half from 
an estimated 2% uncertainty in the relative normalization of 
the FOS and STIS continua. We note that each FOS spectral 
bin is 60 km s _1 wide, so the error in the shift is about half 
an FOS spectral bin. 

The measured change —35 ± 25 km s _1 , equation ©, 
in the position of the shock front between the STIS and 
FOS data is consistent with the predicted change —15.8 ± 
0.7 km s _1 , equation iJTJ. However, the difference is also only 
1.4 standard deviations from no change at all, and it is not 
yet possible to claim a detection. 

3.3 Si III and Si IV features 

As in Paper 1, we have fitted the Si III 1206 A and Si IV 
1394, 1403 A features in the absorption spectrum of SN1006 
to a combination of unshocked and shocked components, the 
two components being assumed to have the same profiles as 
the unshocked and shocked components of the Si II 1260 A 
feature shown in Figure 2] The results are essentially the 
same as found in Paper 1: the Si III and Si IV profiles are 
consistent with being entirely shocked, with no evidence of 
any unshocked component. 

Figure|5|shows, in dotted lines, the fitted Gaussian pro- 
files of shocked Si III and Si IV, the centre and width of the 
Gaussian in each case being fixed equal to those measured 
from the Si II 1260 A feature, as given in Table H The best 
fit ratios of number densities of shocked Si III and Si IV 
relative to shocked Si II are listed in Tabled As elsewhere 
in this paper and in Paper 1, the quoted uncertainties are 
formal la errors, and do not include systematic uncertainties 
associated with the stellar spectrum or with placement of the 
continuum. The systematic uncertainty could well exceed 
the statistical uncertainty. 

If the width of the Gaussian is permitted to be a free 
parameter, the centre being held fixed, then the Si III profile 
prefers a slightly narrower profile, 1070 ± 40 km s -1 , while 
the Si IV profile prefers a slightly broader profile, 1220 ± 
140 km s -1 , but in both cases the width is consistent with 
the 1160 ± 50 km s _1 dispersion of the Si II profile. 

If not only the width but also the centre of the Gaussian 
are permitted to be free parameters (which yields a reliable 
result only for Si IV, since the blue edge of Si III is confused 
by Lya), then the Si IV profile prefers a slightly broader 
width, 1280 ± 140 km s -1 , still consistent with that of Si II, 



but a lower central velocity of 4850 ± 100 km s , which 
is 2-3cr below the 5160 ± 70 km s _1 centre measured from 
Si II, suggesting weakly that Si IV may perhaps extend to 
somewhat lower velocities than Si II. Examination of the 
Si IV line profile, Figure [5] indicates that the fit is pulled to 
lower velocities mainly by an absorption spike at 1405 A (not 
at Si IV 1403 Al). The spike is not so prominent in the FOS 
spectrum, which does however show, like the STIS spectrum, 
a slight excess of Si IV absorption at lower velocities. 

In Paper 1, the Si IV preferred a significantly broader 
width of 1700 ± 100 km s" 1 . The narrower width measured 
here is to a large extent a consequence of dividing by the 
template stellar spectrum of PG 0839+399 - compare the 
raw Si IV profile of the SM star in Figure 0to the corrected 
Si IV profile in Figure |2] 

Since the collisional ionization time of Si II is compara- 
ble to the age of the remnant (Paper 1), it is not surprising 
that collisional ionization of the shocked Si II should lead to 
observable amounts of shocked Si III and Si IV. 

Irradiation of unshocked ejecta by UV and x-ray emis- 
sion from the reverse shock could in principle photoionize 
Si II to Si III and Si IV. Indeed, Paper 1 remarked that the 
ratios of the observed column densities of Si II, III, and IV 
are close to the unshocked ratios pred icted in photoioniza- 
tion tr ials similar to those described bv lHamilton fc SarazirJ 
Jl984h . in which the deflagrated white dwarf model CDTG7 
of Woosley (Woosley 1987, private communicat ion), similar 
to model CDTG5 of IWooslev fc Weaver! (Il987l) . is evolved 
into a uniform ambient medium. Thus the fact that the 
observations point to negligible amounts of unshocked Si III 
and Si IV places constraints on the explosion model - 
for ex ample, the model considered bv lHamilton fc SarazirJ 
( 1984) produces more unshocked Si III and Si IV than 
observed. 

In Paper 1 it was argued that, under a 'simplest' 
set of assumptions, there should be an observable column 
density of blueshifted Si IV, although this prediction is 
sensitive to the assumptions, and is not robust. The STIS 
spectrum, Figure |2] shows no evidence of blueshifted Si IV 
absorption along this sight line, in agreement with previous 
observations. 

3.4 Weaker Si II features 

The absorption spectrum of SN1006 in Figure |2] shows, in 
addition to the broad, redshifted Si II 1260 A feature, weaker 
feat ures of Si II 1304 A and Si II 1527 A. 

iFesen et all Jl988l) first pointed out that the Si II 1527 A 
feature is anomalously weak compared to the principal Si II 
1260 A featu re. The discrepa ncy was confirmed by the FOS 
observations llWu et alJll997f) . and is confirmed again in the 
present observations, Figure [5] 

The high quality of the STIS spectra establishes beyond 
doubt that the discrepancy is real, not just noise. As 
remarked in [J5| the small scale structure in the absorption 
spectrum shown in Figure [5] is real, not noise, and must be 
attributed, at least in unabsorbed (by SN1006) regions of 
the spectra, to differences between the stellar spectra of the 
SM star and the comparison PG 0839+399 star. 

The bottom panel of Figure shows the residual 
spectrum obtained by dividing the absorption spectrum 
by a model in which the weaker Si II 1304 A and Si II 
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1527 A features are assumed to have the same profile as 
the strong Si II 1260 A feature. It is evident that the most 
prominent residual (besides Lya) is around the anomalous 
Si II 1527 A feature. By contrast, the Si II 1304 A feature 
appears to be consistent with Si II 1260 A, in the sense 
that the level of fluctuation in the residual is comparable to 
that in unabsorbed (by SN1006) parts of the spectrum. The 
principal anomaly of the Si II 1527 A feature is that it is 
missing absorption on the red side of the line, by perhaps 
a factor of 2. The Si II 1304 A feature does not share this 
anomaly. 

Note that the red side of the Si II 1304 A feature is 
affected by interstellar C II 1334.5 A, while the blue side 
of the Si II 1527 A feature is affected by interstellar C IV 
1548.2, 1550.8 A. However, these interstellar lines absorb 
only a narrow interval of the high resolution STIS spectra, 
and their presence does not change any conclusions. 

iFesen et all (Il988h proposed that the Si II 1260 A 
feature might be augmented by a contribution from S II 
1250.6, 1253.8, 1259.5 A. Paper 1 argued against this 
proposal on the grounds that that the combined oscillator 
strength of the S II multiplet, 0.005 + 0.011 + 0.016 = 0.033, 
is only l/30th tha t of the Si II 12 60.4 A line, whose oscillator 
strength is 1.007 iMortorJll99lh . This, along with the fact 
that explosion models predict that S and Si are synthesized 
together in a cosmic ratio of order 0.5, suggested that S II 
is likely to contribute no more than 2% to the Si II 1260 A 
feature. 

Could the S:Si ratio in the SN1006 ejecta be substan- 
tially higher than the cosmic ratio of 0.5? Recent studies 
of Type la supernovae have favored off-centre delayed det- 
onation models, in which an explosion begins as a subsonic 
deflagratio m which at som e point develops into a supersonic 
detonation iHoflichl t2006+ Si and S are synthesized in the 
detonation phase by partial burning to nuclear statistical 
equilibrium. The overall S:Si mass ratio is ~ 0.5, but the 
ratio can be higher at lower expansion velocities where 
the burning is more complete. For example, the S:Si ratio 
approaches ~ 2 at an expansion velocity o f 6000 km s _1 
in model CS15DD1 of llwamoto et alJ <ll999T) . Even at such 
enhanced abundance, S would contribute less than 10% to 
the Si II 1260 A feature. 

Since the chief anomaly of the Si II 1527 A feature is 
that it is missing absorption to the red side, an explanation 
in terms of S II absorption would be more plausible if S II 
were mainly unshocked, with little shocked S II. In that case 
S II 1250.6, 1253.8, 1259.5 A would contribute absorption to 
the unshocked component on the red side of the strong Si II 
1260 A feature, but not to the shocked Gaussian component 
on its blue side. However, the collisional ionization rate of 
S II is similar to that of Si II at any electron temperatur e 
exceeding ~ 3 x 10 5 K jLennon et alJll988lfVoronovlll997l) . 
so if S II does contribute to the Si II 1260 A feature, 
then S II should contribute to both unshocked and shocked 
components in about the same amounts as Si II. Thus S II, 
even if present, does not account naturally for the different 
line profiles of the Si II 1527 A and Si II 1260 A features. 
A loophole in the above argument is that if the S:Si ratio 
is strongly varying with radius, from ~ 30 in the unshocked 
gas to < 3 or so in the shocked gas, then it might be possible 
to arrange that S II contributes appreciably only to the 



unshocked component on the red side of the strong Si II 
1260 A feature. But this possibility seems contrived. 

It should be remarked that even if S II does contribute 
appreciably to the Si II 1260 A feature, that would not alter 
the principal conclusion of this paper, that the sharp red 
edge of the strong Si II 1260 A feature indicates a shock 
front. 

Therefore, in the absence of any good alternative, we are 
forced to the conclusion that any discrepancy between the 
Si II lines must arise from differences between the stellar 
spectra of the SM star and the comparison PG 0839+399 
star. Certainly this is a plausible conclusion given the clear 
differences between the spectra at wavelengths unaffected 
by absorption from SN1006 ejecta. It would be useful in 
the future to use theoretical model atmospheres of OB 
subdwarfs to resolve the discrepancies better than we have 
been able to do in this paper. 



4 SUMMARY 

In this paper we have reported a high resolution (5 km s" 1 
FWHM) STIS spectrum of the broad Si II, III, and IV 
features produced b y ejecta in SN1006, seen in abs orption 
against the SM star iSchweizer fc Middleditchlll980ll which 
lies behind and near the projected centre of the remnant of 
SN1006. 

The Si II 1260 A feature, the strongest absorption 
feature, is observed to have an unresolved sharp red edge 
at 7026 + 3 (relative) ± 10 (absolute) km s~\ supporting the 
proposal of Paper 1 that this edge represents the location of 
the reverse shock moving into Si II ejecta. 

The shock model predicts that the position of the 
reverse shock should be changing at the observable rate of 
— 2.7 + 0.1 km s _1 per year. We have attempted to measure 
the change between the position of the reverse shock over 
the 6 year interval between the FOS and STIS data. The 
measured change of —35 ± 25 km s _1 is consistent with the 
predicted change of —15.8+0.7 km s _1 , but is also consistent 
with no change. A future precise detection of the predicted 
change would confirm beyond reasonable doubt that the Si II 
1260 A feature indeed reveals Si II in the process of being 
shocked. 

The pre-shock velocity (7026 km s _1 ), the post-shock 
bulk velocity (5160 km s _1 ), and the post-shock velocity 
dispersion (1160 km s _1 ) of the Si II are all measured from 
the Si II 1260 A profile. Agreement between the measured 
velocities and the shock jump conditions requires that 
almost all the shock energy goes into ions, with little or no 
electron heating. This is consistent with other observational 
evidence that little electron heating occurs in the fast 
shocks in SN1006 ([Lamin g et all Il996l : iGhavamian et al] 
12001 IVink et alJl200SHVinll200^ . 

The STIS spectrum is consistent with the analysis of 
Paper 1, which was based on the lower resolution FOS 
spectrum of the SM star, and the conclusions of that paper 
remain unchanged. 
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